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FLUID CONDITION MONITORING USING 
BROAD SPECTRUM IMPEDANCE SPECTROSCOPY 

Background of the Invention 

[0001] The present invention relates to monitoring the condition of a substance using 
impedance spectroscopy to indicate in real or near real time, i.e., while the substance is being 
used, the physio-chemical condition of a substance based on a correlation of measurements 
from electrical signals, using a statistical technique, to previously determined values. 
[0002] The use of impedance spectroscopy techniques to monitor fluid conditions is 
previously known. For example, U.S. Patent Nos. 6,433,560, 6,380,746, 6,377,052, and 
6,278,281, along with U.S. Published Application 2003/0141882 of Zou et al., all assigned to 
the assignee of the present invention, all teach different configurations of electrodes for 
measuring current and then computing values for impedance. Further, the '281 patent 
teaches comparing the difference in currents at two frequencies with known bulk and 
interfacial impedance measurements to determine fluid conditions. Similarly, the '052 and 
'746 patents teach comparing the difference in impedance values calculated from exciting 
electrodes at two frequencies to determine fluid conditions. 

[0003] The value of using impedance spectroscopy to monitor fluid conditions lies in the 
fact that it is desirable to be able to determine when a fluid, for example, an engine lubricant, 
has degraded to the point where it has either exhausted or come close to exhausting its useful 
life. Similarly, it is desirable to know how many hours of useful life remain with respect to a 
fluid sample. For an application such as monitoring the condition of an engine lubricant, it is 
desirable to be able to monitor the fluid condition while the engine is operating, as opposed 
to performing tests in a laboratory. 

[0004] Lubricating fluids comprise three basic components: (1) base stock, (2) additives, 
and (3) contaminants. These components are known to influence the bulk and interfacial 
properties of the lubricant. Lubricating fluids possess interfacial properties, such as wear 
protection and corrosion protection, which are present at the interface between the fluid and 
the metal it protects. Different portions of the impedance spectrum correspond to bulk and 
interfacial properties. Previous applications of impedance spectroscopy have measured bulk 
properties and interfacial properties separately, but the capability of measuring the two 
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together has not been previously appreciated. Accordingly, it would be desirable if bulk and 
interfacial properties could be measured together in order to give a more complete picture of 
the engine lubricant. 

[0005] Existing methods of monitoring fluid conditions using impedance spectroscopy 
fail to contemporaneously measure a plurality of fluid properties. Rather, prior art methods, 
including those disclosed by the above-mentioned patents and publication, calculate 
impedance in different ways in order to determine a value for one fluid property. It can be 
shown that up to 90% of available information is not utilized when single parameter 
measurement techniques are employed. Thus, there is a need for systems and methods 
capable of implementing a multiple parameter function or functions to make quantitative 
measures of a broad range of fluid condition metrics. 

Brief Summary of the Invention 

[0006] The present invention uses impedance spectroscopy to determine the amounts of 
additives, contaminants, and other components that are present in a fluid. In some 
embodiments, the present invention is applied to engine lubricants. 

[0007] The invention comprises using the impedance spectrum to monitor the condition 
of a fluid in real time by disposing an electrode mechanism in the fluid and exciting the 
electrode mechanism sequentially with a specified number of alternating voltages, wherein 
each of the alternating voltages is at a different frequency in a range of frequencies. The 
invention measures the current in the electrode mechanism at each of said frequencies in the 
range of frequencies. From the measurement of current at each of the frequencies, resistive 
and reactive impedance values at each of the frequencies are calculated and used to predict at 
least one property in the fluid in situ. 

[0008] In one embodiment, the invention comprises a method of analyzing a substance, 
comprising the steps of: (1) generating first and second plots of spectra over a range of 
frequencies; (2) creating a third plot that comprises the first and second plots by sequentially 
assigning x-values to selected frequencies in the first plot and selected frequencies the second 
plot; and (3) repeating steps (1) - (2) at least once to generate a plurality of third plots. Next, 
the method comprises building a spectral matrix that comprises data taken from the plurality 
of third plots and a result matrix comprising known quantities of a plurality of components in 
the substance. A Principal Component Analysis is performed on the spectral matrix to 
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identify at least one principal component having significant influence on the spectral matrix. 
A reduced spectral matrix having at least one column, wherein each column in the reduced 
spectral matrix is associated with a principal component having significant influence on the 
spectral matrix, is next created. Then, a statistical technique is performed that uses the 
reduced spectral matrix and the result matrix to create at least one prediction equation. The 
at least one prediction equation is then used to predict at least one property in a second 
substance in situ. 

[0009] In some embodiments, at least one datum derived from a Nyquist plot is added to 
the spectral matrix. The at least one datum derived from a Nyquist plot may include at least 
one datum from the bulk region of the Nyquist plot and at least one datum from the 
interfacial region of the Nyquist plot. 

Brief Description of the Drawings 

[0010] FIG. 1 depicts a system by which predictions of fluid conditions are made. 
[0011] FIG. 2 describes a system that may be used for collecting data used to develop an 
information library used to predict fluid conditions. 

[0012] FIG. 3 depicts a flowchart providing an overview of the method by which an 
information library is developed and used to predict fluid conditions in real time. 
[0013] FIG. 4 provides a flowchart describing a method for building a matrix of spectral 
data, i.e., a spectral matrix. 

[0014] FIG. 5 provides a flowchart describing a method for building a result matrix. 
[0015] FIG. 6 shows the process by which the present invention uses Principal 
Component analysis (PCA). 

[0016] FIG. 7 shows Bode plots of resistive and reactive impedance spectra 
superimposed on the same graph. 

[0017] FIG. 8 shows an example of combined Bode plots on a graph. 

[0018] FIG. 9 provides an example of a result matrix populated with data from laboratory 

tests representing conditions of an engine lubricant. 

[0019] FIG. 10 shows regression coefficients plotted on a graph. 

[0020] FIG. 1 1 shows a reduced spectral matrix. 

[0021] FIG. 12 shows a chart giving examples of data, derived from a Nyquist plot of 
resistive impedance versus reactive impedance, that can be added to a combined Bode plot. 
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[0022] FIG. 13 shows a graph comprising three Nyquist plots for three different 
lubricating fluid samples. 

[0023] FIG. 14 describes an alternative embodiment in which data from each of the Bode 
plots of resistive and reactive impedance spectra are processed separately, with results then 
combined for use in predicting fluid properties. 

Detailed Description of the Invention 

[0024] SYSTEM OVERVIEW 

[0025] Figure 1 provides an overview of an exemplary system 100 by which the 
condition of a substance such as a fluid is determined in real or near real time. Although the 
invention is described herein with reference to system 100, those skilled in the art will 
appreciate that other configurations and other components could support the claimed systems 
and methods so long as impedance spectroscopy is used to generate values for resistive and 
reactive impedance, which values could be used to perform calculations as described below 
to predict fluid conditions in situ and generally while the fluid is being actively used. 
[0026] Electrode mechanism 102 is disposed in fluid 104 contained in vessel 106. The 
electrode mechanism may be any of those that are known in the art of impedance 
spectroscopy, including those described in aforementioned U.S. Patent Nos. 6,433,560, 
6,380,746, 6,377,052, and 6,278,281, and U.S. Published Application 2003/0141882 of Zou 
et al., all of which are hereby incorporated herein by reference. Present practice is to use the 
cylindrical probe mechanism disclosed in U.S. Published Application 2003/0141882 of Zou 
et al. In some embodiments of the present invention fluid 104 is an engine lubricant and 
vessel 106 is an engine crankcase. 

[0027] Power supply 1 14 powers oscillator 1 12, a computing system or device such as 
microcomputer 116, and excitation driver 108. Excitation driver 108 receives an input from 
an oscillator 112. In one embodiment, oscillator 112 provides voltage in a range from 
approximately 75 kilohertz and 0.0075 hertz. However, those skilled in the art will 
appreciate that the invention could be practiced using voltages at frequencies higher than 75 
kilohertz and/or lower than 0.0075 hertz. In one embodiment, oscillator 1 12 sequentially 
provides voltages at a specified number of different frequencies to excitation driver 108. 
Excitation driver 108, upon receiving input from oscillator 112, excites electrode mechanism 
102, and current sensor 1 18 provides input to microcomputer 1 16. Present practice is for 
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excitation driver 108 to provide voltages of approximately 250 millivolts RMS, although 
those skilled in the art will appreciate that the invention could be practiced using voltages at 
other levels. Temperature sensor 1 10 measures the temperature of fluid 104, and provides 
input to microcomputer 116. 

[0028] Microcomputer 116 calculates and stores in memory values for resistive and 
reactive impedance, not shown in Figure 1, based on the input from current sensor 118. 
Resistive and reactive impedance are sometimes referred to as real and imaginary impedance, 
respectively. Fluid temperature from temperature sensor 1 10 is input to computer 1 16. 
Baseline values 122, comprising previously predicted values of fluid properties, are stored in 
information library 120. Each of the baseline values 122 has as attributes the identity of the 
one fluid property with which it is associated, a fluid temperature, and a value representing 
the expended useful life of the fluid. Fluid temperature is usually expressed in degrees 
Celsius, and the value representing expended useful life may be expressed in hours. Baseline 
values 122 are determined in a laboratory external to system 100, and comprise expected 
fluid property values for a fluid of a given age at a given temperature. In some embodiments 
baseline values 122 comprise expected fluid property values for a fluid that has been used in 
an engine for a given number of hours, and is currently at a given temperature. 
[0029] Microcomputer 1 16 uses information library 120 to determine the condition of 
fluid 104 by using the temperature input from temperature sensor 110 and the calculated 
values for resistive and reactive impedance to retrieve the correct prediction equations 124 
from information library 120. Prediction equations 124 are described in detail below. 
Microcomputer 116 then uses prediction equations 124 to calculate fluid property values 128. 
Microcomputer 116 then obtains one or more property predictions 126 of values relating to 
the condition of fluid 104 by comparing at least one predicted fluid property value 128 to at 
least one of baseline values 122. Property predictions 126 will be stored in the memory of 
microcomputer 116, and may optionally be stored in information library 120, although this 
optional configuration is not depicted in Figure 1. 

[0030] In some embodiments, information library 120 is a component of microcomputer 
1 16. In some embodiments, microcomputer 1 16 is further capable of providing as output an 
end of life (EOL) measurement, a remaining useful life (RUL) measurement, or both. In 
these embodiments, baseline values 122 are used in conjunction with the results from using 
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the prediction equations 124, that is, fluid property values 128, to make a determination of 
whether fluid 104 is at or near the end of its useful life. Methods by which EOL and RUL 
determinations are made are discussed in more detail below. 

[0031] Figure 2 describes a system 200 that may be used for collecting data used to 
develop information library 120, including baseline values 122. Although the invention is 
described herein with reference to system 200, those skilled in the art will appreciate that 
other configurations and other components could support the claimed systems and methods 
so long as impedance spectroscopy is used to generate values for resistive and reactive 
impedance. System 200 includes a sample fluid 204 at a predetermined, constant 
temperature contained in a vessel 206. Electrode mechanism 202 is disposed in sample fluid 
204. In one embodiment, sample fluid 204 is an engine lubricant. 

[0032] Power supply 214 powers oscillator 212, a computing system such as computer 
216, and excitation driver 208. Excitation driver 208 receives a sequence of inputs from 
oscillator 212. Oscillator 212 provides voltage at a specified number of frequencies in a 
range of frequencies. As currently practiced the invention uses frequencies in a range from 
approximately 75 kilohertz and 0.0075 hertz. However, those skilled in the art will 
appreciate that the invention could be practiced using voltages at frequencies higher than 75 
kilohertz and/or lower than 0.0075 hertz. In one embodiment, the specified number of 
frequencies is 71. For each input from oscillator 212, excitation driver 208 excites electrode 
mechanism 202, and current sensor 218 provides input to computer 216. Present practice is 
for excitation driver 208 to provide voltages of approximately 250 millivolts RMS, although 
those skilled in the art will appreciate that the invention could be practiced using voltages at 
other levels. Temperature sensor 210 measures the temperature of fluid 204, and provides 
input to computer 216. Computer 216 generates calculated values 228 for resistive and 
reactive impedance based on the input from current sensor 218. Computer 216 stores the 
calculated values 228 in database 220 such that calculated values 228 are associated with the 
specified temperature and the frequency in the range of frequencies for which the calculated 
values for resistive and reactive impedance were obtained. In some embodiments, database 
220 is a component of computer 216. 

[0033] Computer 216 may be further configured to generate prediction equations 224 
according to the method described in detail below. As part of this method, baseline values 
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222 are used to populate a result matrix, as described below. Prediction equations 224 may 
be stored in database 220. At some point prior to use of system 100 for prediction of fluid 
conditions in real time, prediction equations 224 may be copied from database 220 into 
information library 120 in system 100, wherein prediction equations 124 then comprise 
prediction equations 224. Similarly, baseline values 222 may be copied to baseline values 
122. 

[0034] OVERVIEW OF A METHOD FOR PREDICTING FLUID PROPERTIES 
[0035] Figure 3 depicts a flowchart providing an overview of the method by which 
information library 120, including prediction equations 124 and baseline values 122, is 
developed and used to predict fluid conditions in real time. When reviewing Figure 3 it is 
helpful to bear in mind that the overall goal of the process being described is to develop a set 
of prediction equations 124, each of which can be used to predict a fluid property. These 
prediction equations 124 will take the form: 

Y u = b 0 + bjXj + b 2 X 2 + . . . + b f X f (1) 
[0036] where Yi j represents the j th fluid property out of a given number of fluid 
properties being measured in the i th sample out of n fluid samples being used, each Xi . . . Xf 
represents a value associated with an impedance reading at a particular frequency, and each 
bo ... bf is what is called a loading coefficient. One goal of the inventive method is to solve 
for the loading coefficients using values for Xi . . . Xf determined from a process using 
impedance spectroscopy, certain statistical techniques, and baseline values 122 for Yjj that 
are known from prior laboratory testing. The developed equations, i.e., prediction equations 
124, can then be used to determine a set of values Yy in real time, i.e., the property 
predictions 126 described above with reference to Figure 1. 

[0037] Turning now to Figure 3, block 300 represents the process of assembling a matrix 
of spectral data relating to a set of n fluids. In one embodiment, each of the n fluids will 
represent fluid conditions for the fluid at n different values for expended useful life. The 
goal of the process represented by block 300 is to develop a matrix of data, hereinafter 
referred to as the spectral matrix, representing impedance spectra obtained using each of the 
n fluids that can then be subjected to statistical analysis. The form of the spectral matrix is as 
follows: 
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^0.0^1,0 ^0,1^1.1 ^O.p.^O.p-l 

^1.0^1.0 ^1.1^1,1 

The dimensions of the spectral matrix will be n rows by p columns, where n is the number of 
fluid samples, as discussed above, and p is related to, and sometimes equal to, the number of 
frequencies for which calculations of impedance have been used. In one embodiment p is 
equal to one-hundred and forty-two, and represents two times the number of frequencies for 
which impedance calculations have been used. The process of creating the spectral matrix is 
described in detail below with reference to Figure 4. 

[0038] Block 302 represents the process of assembling a result matrix for the n fluids, 
i.e., the same set of n fluids for which the spectral matrix was constructed in Block 300 
above. As noted above baseline values 222 comprise data used in the result matrix. The data 
in the result matrix is achieved from traditional analytical laboratory procedures for 
measuring fluid properties. While by no means limited to the following, examples of the 
laboratory tests used to determine values for the result matrix include ASTM D-445 (40 or 
100 degree vis), ASTM D-4739 (TBN), ASTM D-2869 (TBN, i.e., Total Base Number), 
ASTM D-664(TAN, i.e., Total Acid Number), ASTM D_5967 (per cent SOOT), and ASTM 
D-5185 (ICP Elementals). The form of the result matrix is as follows: 

Y Y 

Y Y 

Y Y 

The data in the result matrix represent fluid conditions with respect to f fluid properties in n 
sample fluids. Accordingly, the dimensions of the result matrix are n rows by f columns. 
The result matrix is used to find the loading coefficients for the prediction equations 124 (1), 
as discussed in more detail below. The process of creating the result matrix is described in 
further detail below with reference to Figure 5. 

[0039] Returning to Figure 3, block 304 represents the optional process of applying pre- 
processing functions to the spectral matrix before subjecting the spectral matrix to statistical 
analysis. The pre-processing functions used in the present invention are all well known, and 

include but are not limited to mean-centering, taking the first or second derivative of the data, 
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smoothing the data, sample averaging, and differencing. Differencing comprises taking the 
difference between two values related to a fluid sample, such as the difference between 
values for resistive and reactive impedance, or the difference between data values for a new 
fluid and a fluid that has been used, e.g., a new engine lubricant and a lubricant that has been 
used in an engine operating for a number of hours. 

[0040] Block 306 represents the process of performing a Principal Component Analysis 
(PCA) on the spectral matrix. PCA is a technique for analyzing a set of data to determine 
underlying independent factors that influence the data. Applied in the context of the present 
invention, PCA provides the advantage of reducing the number of variables in the spectral 
matrix to a set of variables for which there are as few common variations in the data as 
possible. By using PCA to create a set of principal components that represent the major 
changes in the impedance spectra found in the spectral matrix, the present invention creates a 
simplified spectra that can be subjected to meaningful statistical analysis that would not be 
practical with the entire spectral matrix. The advantages of PCA can be seen by noting that 
in current practice of the present invention the number of columns in the spectral matrix is 
often reduced from one-hundred and forty-two to less than ten, allowing the creation of a set 
of meaningful data, wherein the size of the data set is such that it is practical to apply 
statistical techniques to the data. PCA in the context of the present invention is discussed in 
more detail below with reference to Figure 6. 

[0041] Block 308 represents the process of creating a reduced spectral matrix by 
redeveloping the spectral matrix with principal components selected from the output of the 
PCA. Specifically, as discussed below with reference to Figures 6 and 10, principal 
components are selected for the reduced spectral matrix if they show significant influence on 
the prediction of fluid properties. The selected principal components are then associated with 
the X values from the spectral matrix that correspond to them, and placed in the reduced 
spectral matrix. The process of creating a reduced spectral matrix is described in detail 
below with reference to Figure 6. 

[0042] Block 310 represents the process of performing a statistical analysis, using the 
result matrix and the reduced spectral matrix, to find loading coefficients that can be used in 
the prediction equations 124 (1) described above. Some embodiments of the present 
invention use Principal Component Regression (PCR) to determine the loading coefficients, 
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while other embodiments use Multivariate Least Squares Regression (MLR) to determine the 
loading coefficients. PCR and MLR are both well known. Still other embodiments use well 
known nonlinear regression methods such as Group Methods of Data Handling. 
Additionally, other known methods of statistical analysis could be used to determine the 
loading coefficients. The use of statistical techniques in the context of the present invention 
is discussed in more detail below. 

[0043] Block 312 represents the process of using the prediction equations 124 (1) in 
system 100 of Figure 1 to obtain predicted property values 128. Predicted property values 
128 are in turn compared with baseline values 122 to obtain property predictions 126. This 
processing may be done by microcomputer 116. In one embodiment, predicting fluid 
properties for a fluid in situ comprises determining the likely expended useful life of an 
engine lubricant by determining whether property predictions 126 either exceed or fail to 
meet predetermined threshold values, which determination may be done by microcomputer 
1 16. From this determination an estimate of Remaining Useful Life (RUL) of the lubricant 
can be made and output using microcomputer 1 16. Similarly, if it is the case that, based on 
the determination of the likely age of the fluid, the fluid is near the end of its useful life, an 
End of Life (EOL) determination can be made and output from microcomputer 1 16. The 
procedure for making an RUL or EOL estimation in the context of the present invention is 
discussed in more detail below. 
[0044] CREATION OF A SPECTRAL MATRIX 

[0045] Figure 4 provides a flowchart describing an approach for building a matrix of 
spectral data. This approach is described herein with reference to system 200, although it is 
to be appreciated that other systems for conducting impedance spectroscopy, and recording 
and performing measurements and calculations therefrom, could be used. In step 400, a 
sample fluid 204 having a known expended useful life is provided at a specified temperature. 
In step 402, the expended useful life in hours and the specified temperature are recorded in 
database 220. It will be appreciated that recording the specified temperature is important 
inasmuch as many fluid measurements are a function at least in part of temperature. 
Similarly, many fluid measurements are a function at least in part of expended useful life. 
One of the objects of the present invention is to be able to determine a value for the expended 
useful life of a fluid when temperature is known. 
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[0046] In step 404, with reference to system 200, calculated values for resistive and 
reactive impedance 224 associated with predetermined properties of sample fluid 204 are 
recorded in database 220. Baseline values 222 are recorded in database 220 such that they 
are associated with the specified temperature of sample fluid 204 as well as the known age in 
hours of sample fluid 204. The predetermined fluid properties associated with baseline 
values 222 may include bulk and/or interfacial properties of the sample fluid 204 comprising 
amounts of additives such as zinc dithiodialkylphosphates (ZDDPs) and contaminants such 
as soot as well as interfacial properties such as wear protection. Baseline values 222 will 
have been previously determined through one of a variety of testing means known to those 
skilled in the art, as discussed above. 

[0047] In step 406, test parameters are selected, comprising a range of frequencies to be 
tested, a specified number of specific frequencies in the range of frequencies to be tested, and 
the specific frequencies that will be tested. Present practice is to select 71 frequencies in the 
range from 75 kilohertz to 0.0075 hertz. A listing, in hertz, of frequencies used in at least 
one instance of present practice of the invention is as follows: 75000, 59574.62, 47321.8, 
37589.04, 29858.04, 23717.08, 18839.15, 14964.47,11886.7, 9441.941, 7500, 5957.462, 
4732.18, 3758.904, 2985.804, 2371.708, 1883.915, 1496.447, 1188.67, 944.1941, 750, 
595.7462, 473.218, 375.8904, 298.5804, 237.1708, 188.3915, 149.6447, 118.867, 94.41941, 
75, 59.57462, 47.3218, 37.58904, 29.85804, 23.71708, 18.83915, 14.96447, 11.88670, 
9.44194, 7.50000, 5.95746, 4.73218, 3.75890, 2.98580, 2.37171, 1.88391, 1.49645, 1.18867, 
0.94419, 0.75000, 0.59575, 0.47322, 0.37589, 0.29858, 0.23717, 0.18839, 0.14964, 0.11887, 
0.09442, 0.07500, 0.05957, 0.04732, 0.03759, 0.02986, 0.02372, 0.01884, 0.01496, 0.01189, 
0.00944, and 0.0075. It should be understood that the invention is not limited to any 
particular frequencies, number of frequencies, or range of frequencies, and that the above list 
of frequencies is given for illustrative purposes only. The above-listed frequencies were 
chosen because they are somewhat evenly distributed across the impedance spectrum and 
have been found to yield good laboratory results. 

[0048] In step 408, again referring back to system 200, oscillator 220 inputs a sequence 
of voltages at each of the specific frequencies selected to excitation driver 208, which sends 
current at each of the specific frequencies to electrode mechanism 202. Current sensor 218 
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detects the amount of current at each of the specific frequencies from electrode mechanism 
202, and sends an amount of current as an input to computer 216. 

[0049] In step 410, computer 216 calculates resistive and reactive impedance values for 
each of the specific frequencies based on the current at each of the specific frequencies, and 
stores the calculated values 224 in database 220. 

[0050] In step 412, a Bode plot of the resistive impedance spectra is created. Bode plots 
are well known in the art. A Bode plot comprises a two-dimensional graph in which the x 
axis is comprised of the logarithms of the frequencies against which impedance is plotted, 
and the y axis is comprised of values for impedance. In step 414, a Bode plot of the reactive 
impedance spectra is created. 

[0051] Figure 7 shows Bode plots of resistive and reactive impedance spectra 
superimposed on the same graph 700. X axis 702 comprises values for the logarithms of 
frequencies in the impedance spectra. Y axis 704 comprises calculated values for impedance 
corresponding to the frequencies whose logarithms are represented in x axis 702. Plots 706 
comprise plots of the resistive impedance spectrum for different fluid samples. Plots 708 
comprise plots of the reactive impedance spectrum for different fluid samples. 
[0052] Returning to Figure 4, in step 416 the values on the x axis of the Bode plot of the 
resistive impedance spectra are replaced with integers ranging from zero to a number one less 
than the specified number of frequencies that were selected in step 406. In step 418, the 
values on the x axis of the Bode plot of the reactive impedance spectra are replaced with 
integers ranging from the specified number of frequencies that were selected in step 406 to a 
number that is one less than two times the specified number of frequencies. For example, if 
the specified number of frequencies was seventy-one, values on the x-axis would range from 
zero to one-hundred and forty-one. As discussed below, the number of columns p in the 
spectral matrix is the number identified in the present step that is two times the specified 
number of frequencies. Present practice is to include one-hundred and forty-two columns in 
the spectral matrix. 

[0053] In step 420, the Bode plot of the resistive impedance spectra created in step 412 
and modified in step 416 is combined with the Bode plot of the reactive impedance spectra 
created in step 414 and modified in step 418. In one embodiment, the Bode plots of the 
resistive and reactive spectra are combined so that the plots are laid "head to tail", with the 
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maximum plotted value on the x axis of the plot of the resistive impedance spectra one unit 
to the left of the minimum value on the x axis of the plot of the reactive impedance spectra. 
It should be understood that the order in which the Bode plots of resistive and reactive 
impedance are combined is not essential to the practice of the invention. For convenience, 
the approach and examples discussed herein place a Bode plot of resistive impedance to the 
left of a Bode plot of reactive impedance on a combined graph. This order could just as 
easily be reversed. In fact, points from the Bode plots of the resistive and reactive spectra 
could be placed on the combined plot in any order without making a difference to the results 
achieved by practicing the invention. 

[0054] Figure 8 shows an example of combined Bode plots on a graph 800. Y axis 804, 
which comprises the same scale as found on y axis 704, comprises calculated values for 
impedance corresponding to the frequencies whose logarithms are represented in x axis 702. 
With respect to x axis 802, however, the logarithms of the frequencies in the impedance 
spectra have been replaced by integers ranging from zero to one-hundred and forty-one. 
Resistive impedance spectra 806 lie above the x axis values less than or equal to seventy. 
Reactive impedance spectra 808 lie above the x axis values greater than or equal to seventy- 
one. 

[0055] Returning to Figure 4, in step 422 a test is done to determine whether there are 
any other data points to be added to the combined plot. It has been found that adding certain 
data points to the combined plot can increase the ability of the invention to predict fluid 
properties. Examples of data points that may be added to the combined plot are discussed 
below with respect to Figures 12 and 13. 

[0056] If there are other data points to be added to the combined plot, then they are added 
in step 424 by placing them either to the left of the minimum x value presently graphed on 
the combined plot, or to the right of the maximum x value currently graphed on the combined 
plot. If there are not any data points to be added to the combined plot, then the process 
continues to step 426. 

[0057] In step 426, a test is done to determine if there is another sample fluid 204 to be 
tested. If the answer is yes, the method returns to step 400. If the answer is no, the method 
proceeds to step 428. In general, there will be n iterations of the method, n being the number 
of fluid samples at different ages that are being tested. 
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[0058] In step 428 a spectral matrix, taking the form described above, is constructed. As 
discussed above, the spectral matrix will have n rows and p columns. The number p 
generally represents the number of frequencies for which impedance readings are included in 
the matrix. However, as mentioned above, in one embodiment of the present invention, p is 
actually twice the number of frequencies for which impedance measurements were taken 
because resistive impedance and reactive impedance are placed together on combined plots. 
For example, in one embodiment impedance measurements are taken at seventy-one different 
frequencies, and p is therefore one-hundred and forty-two when resistive and reactive 
impedance plots are combined. Each number X, or x value, in the spectral matrix represents 
an integer from zero to p, and each number Y represents an impedance value on the 
combined graph associated with its corresponding x value. 

[0059] ADDING DATA TO THE COMBINED PLOT AND SPECTRAL MATRIX 
[0060] Figure 13 shows graph 1300 comprising three Nyquist plots of the impedance 
spectra for three different lubricating fluid samples. X-axis 1302 comprises values for 
resistive impedance, also known as real impedance, denoted Z\ Y-axis 1304 comprises 
values for reactive impedance, also known as imaginary impedance, denoted Z". As will be 
understood by those skilled in the art, for each of the three plots on graph 1300 the data 
points to the left of the minimum value for Z M represent the fluid bulk. Likewise, those 
skilled in the art will understand that the data points to the right of the minimum value for Z M 
for each of the three plots on graph 1300 represent the interfacial region between the fluid 
and metal. In one embodiment, the interfacial region would represent where the fluid is in 
contact with an engine. Some of the data points discussed below are not contained on the 
graph 1300 but will be known from data recorded from generating the impedance spectra. 
[0061] Figure 12 depicts chart 1200 giving examples of data that can be added to the 
combined plot derived from a Nyquist plot of resistive impedance versus reactive impedance. 
Column 1202 lists the identifiers for each of the samples for which data is provided. Note 
that the data in chart 1200 was taken from the same fluid samples as was used to generate the 
Bode plots shown in Figures 7 and 8. Further, although graph 1300 shows only three 
Nyquist plots, in actual practice of the invention a Nyquist plot would be generated for each 
of the data samples in column 1202. 
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[0062] Column 1204 represents the Z' value where Z" is minimum. Column 1206 
represents the Z M value where Z M is minimum. As noted above, the point at which Z" is 
minimum denotes the boundary between the bulk and interface regions of the Nyquist 
spectrum. 

[0063] Column 1208 represents the frequency in the spectra at which Z" is minimum. 
Graph 1300 does not represent this value. However, all of the values for Z M will be contained 
in the spectral matrix, and moreover, it is possible to determine from the spectral matrix the 
frequency in the impedance spectrum at which each value was recorded, inasmuch as the 
Bode plots used to create the spectral matrix originally contained the logarithm of frequency 
on the x axis. 

[0064] Column 1210 represents the maximum Z' value within the total data set. In some 
cases this value is the "Nyquist Max," that is, the Z ! value between the bulk and interfacial 
regions of the Nyquist spectrum. In other cases this value could be the Z' value for data 
points associated with the lowest or highest frequencies. 

[0065] Column 1212 represents the minimum Z' value within the total data set. In some 
cases this value is the "Nyquist Min," that is, the Z' value between the bulk and interfacial 
regions of the Nyquist spectrum. In other cases this value could be the T value for data 
points associated with the lowest or highest frequencies. 

[0066] Columns 1214, 1216, 1218, and 1220 all contain data associated with the points in 
the bulk region of the Nyquist spectrum, that is, the points between the origin of the graph 
and the minimum value for Z". As is known, these points describe a semicircle. Certain 
information about the semicircle, or the circle that would result from completing the 
semicircle, can be helpful in predicting fluid properties. 

[0067] Column 1214 represents the Z' value for the centerpoint of the circle in the bulk 
region of the Nyquist spectrum, that is, the circle completed by the semicircle drawn from the 
leftmost point on the x axis on which data is plotted to the point on the x axis at which Z" is 
plotted. This circle is sometimes referred to as the bulk circle. 

[0068] Column 1216 represents the Z" value for the centerpoint of the centerpoint of the 
bulk circle. 
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[0069] Column 1218 represents a measurement in radians of the angle between the x axis 

and a line drawn through the origin of the graph and the centerpoint of the bulk circle. This 

measurement is referred to as the depression angle of the bulk circle. 

[0070] Column 1220 represents a calculation of the radius of the bulk circle. 

[0071] Columns 1222, 1224, 1226, and 1228 all contain data associated with the points in 

the interfacial region of the Nyquist spectrum, that is by the points to the right on the x axis 

of the minimum value for Z". As is known, these points describe a semicircle. Certain 

information about the semicircle, or the circle that would result from completing the 

semicircle, can be helpful in predicting fluid properties. 

[0072] Column 1222 represents the Z' value for the centerpoint of the circle in the 
interfacial region of the Nyquist spectrum, that is, the circle completed by the semicircle 
drawn from the points to the right on the x axis of the minimum value for Z". This circle is 
sometimes referred to as the interface circle. 

[0073] Column 1224 represents the Z" value for the centerpoint of the interface circle. 
[0074] Column 1226 represents a measurement in radians of the angle between the x axis 
and a line drawn through the origin of the graph and the centerpoint of the interface circle. 
This measurement is referred to as the depression angle of the interface circle. 
[0075] Column 1228 represents a calculation of the radius of the interface circle. 
[0076] The data in Figure 12 can be added to the spectral matrix and used to predict fluid 
properties either by itself or alongside other values, such as the selected values for resistive 
and reactive impedance already discussed. This data can be included in the spectral matrix 
by placing it on a graph either to the left or to the right of the combined Bode plots discussed 
above. Some, only one, or all of the data points discussed above may be used. In addition, 
skilled artisans will recognize that adding other data related to a scan of the impedance 
spectrum is within the scope and spirit of the present invention. 
[0077] CREATION OF THE RESULT MATRIX 

[0078] Figure 5 provides a flowchart describing a method for building a result matrix. In 
step 500, a fluid sample having a known expended useful life is tested in a laboratory to 
determine values for a predetermined set of fluid conditions, and these values are recorded. 
As discussed above with respect to Figure 2, these values may be recorded in database 220 as 
baseline values 222. Some of the laboratory tests used in present practice are mentioned 
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above with reference to block 302 in Figure 3. In step 502, a check is done to determine 
whether there are more fluid samples to be tested. If the answer is yes, the process returns to 
step 500. If the answer is no, the method proceeds to step 504. In general, there will be n 
iterations of the method, n being the number of fluid samples having different expended 
useful lives that are being tested. 

[0079] In step 504 a result matrix, taking the form described above with reference to 
block 302 of Figure 3, is constructed. Each number Yij in the result matrix represents a fluid 
condition with respect to the jth fluid property from a total of f fluid properties tested in the 
ith fluid sample from a total of n fluid samples. 

[0080] Figure 9 provides an example of a result matrix 900 populated with empirically 
achieved data representing properties of an engine lubricant. Sample column 902 contains an 
identifier for each of the fluids represented in the result matrix. Hours column 904 contains 
the number of hours during which the lubricating fluid has been in operation in the engine, 
i.e., the expended useful life of the lubricating fluid. It should be understood that all of the 
data in the results matrix represents fluid conditions at a single temperature. Note that each 
value in column 904 contains the suffix "_70", which indicates that the data in result matrix 
900 represents fluid conditions at seventy degrees Celsius. 

[0081] Result columns 906, 908, 910, 912, 914, 916, and 918 each contain values 
representing the measurement of a particular fluid property determined by performing 
laboratory tests as described above. Result matrix 900 depicts fifteen sample lubricants in 
which seven fluid properties have been measured. Accordingly, when put into the form 
given above, result matrix 900 is a fifteen by seven matrix; that is, with respect to result 
matrix 900, n is fifteen, because 15 sample fluids have been analyzed, and f is seven. 
[0082] PRINCIPAL COMPONENT ANALYSIS 

[0083] Figure 6 shows the process by which the present invention uses Principal 
Component analysis (PCA). PCA is a well known technique, described in a plethora of 
literature, including the articles "Principal Component Analysis Methods" and "Discriminant 
Analysis, The PCA/MDR Method", both found on the website of Thermo Galactic of Salem, 
New Hampshire at "www.galactic.com/Algorithms/pca.htm" and 

"www.galactic.com/Algorithms/discrim_pca.htm" respectively and incorporated by reference 
herein. PCA is also explained in the following, all of which are incorporated by reference 
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herein: Michael Palmer, "Principal Component Analysis" found on the world wide web at 
"www.okstate.edu/artsci/botany/ordinate/PCA.htm"; StatsSoft, Inc., "Principal Components 
and Factor Analysis" found on the world wide web at 

"www.statsoftinc.com/textbook/stfacan.html"; "Principal Component Analysis" found on the 
world wide web at "www.casaxps.cwc.net/FactorAnalysis.htm". Because PCA is well 
known, it will be described herein only to the extent necessary to explain how PCA is applied 
in the context of the present invention. 

[0084] The goal of PCA is to reduce the number of elements in a data set by selecting 
principal components that are associated with the most variation in the data set. PCA 
comprises iteratively removing independent variations from a data set. Thus, in the context 
of the present invention the goal of PCA is to create a reduced spectral matrix that comprises 
a subset of the spectral matrix discussed above with respect to Figure 4. PCA operating on 
spectral data uses the following relationship: 

A = SF + E A (2) 
where A is the n by p spectral matrix, S is an n by f score matrix, F is the f by p matrix 
containing the principal components, i.e., a principal components matrix, and E A is an n by p 
error matrix. F is sometimes referred to as a matrix of eigenvectors because it is used to 
recreate the spectral matrix. As above, n is the number of spectral samples in the spectral 
matrix and p is the number of data points represented in the spectra. The number f here 
represents the number of principal components. Thus, PCA applied to spectral data depends 
on the theory that the expression SF can be used to recreate spectral matrix A. 
[0085] The actual calculations required for PCA in the present invention are performed 
by a statistical software package. Software packages that have been used to practice the 
invention include The Unscrambler® from Camo Technologies of Woodbridge, New Jersey; 
Spectrum Quant+™ from PerkinElmer, Inc. of Wellesley, Massachusetts; and MatLab® 
from Mathworks, Inc. of Natick, Massachusetts. Those skilled in the art will appreciate and 
understand the usage of such packages. Accordingly, the process described in Figure 6 is 
relatively simple. At step, 600, the spectral matrix is provided as an input to the PCA 
process. At step 602, the aforementioned statistical software package provides output of the 
PCA process in the form of a principal component matrix. 
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[0086] At step 604, a set of regression coefficients associated with the principal 
components matrix is generated. This step serves the purpose of enabling selection of the 
principal components that appear to have a significant influence on the spectral matrix once 
the principal components are identified in the principal component matrix. That is, it is 
necessary to determine which principal components will be useful in recreating the spectral 
matrix, i.e., which principal components have significant influence over variations in the data 
set. 

[0087] There are different known methods of identifying principal components, and the 
particular method used to select principal components is not critical to practicing the 
invention. However, present practice is to select principal components by generating a set of 
regression coefficients and plotting the regression coefficients over the points in the 
impedance spectra; this plot is called a regression spectrum. Regression coefficients are 
simply the set of coefficients obtained by regressing values in the impedance spectra against 
values representing known fluid properties, i.e., selected values from the result matrix. One 
skilled in the art will appreciate that standard statistical software packages, such as those 
discussed above, and a variety of regression methodologies, could be used to generate the 
regression spectrum. 

[0088] Step 604 can be made clear by way of example. One fluid property measured in 
engine lubricants is Total Base Number (TBN). A predicted value for TBN with respect to a 
fluid sample at a given temperature and age can be retrieved from the result matrix, discussed 
above with reference to Figures 5 and 9, and regressed against the original impedance spectra 
generated as part of the process discussed above with reference to Figure 4. Performing this 
regression generates a set of regression coefficients for TBN. These regression coefficients 
can then be plotted as shown in plot 1006 on graph 1000 of Figure 10. X axis 1002 
comprises integers in the range of the integers contained on the x axis of the combined Bode 
plot discussed above with reference to Figures 4 and 8. Consistent with one embodiment 
referenced in examples previously given, plot 1006 is plotted from points zero to one- 
hundred and forty-one in relation to x axis 1002. Y axis 1004 comprises a range of values of 
coefficients generated from the regression of TBN against the impedance spectra. 
[0089] A visual inspection of plot 1006 reveals that positive or negative peaks occur at 
points 18, 25, 52, 63, 70, 85, 92, 129, and 141 lying above x axis 1002. That is, the 
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regression coefficients at these points have a significant magnitude, and these regression 
coefficients thus exert a relatively higher degree of influence over the prediction of a value 
for TBN than other regression coefficients in plot 1006. Accordingly, columns 
corresponding to points 18, 25, 52, 63, 70, 85, 92, 129, and 141 are selected from the spectral 
matrix and placed in the reduced spectral matrix in step 606. 

[0090] It should be noted that, with respect to the above example, the present invention 
does not require using the regression spectrum for TBN, as opposed to other potentially 
available regression spectra, to pick points to be included in the spectral matrix. Other fluid 
properties, either individually as was the case with TBN in the above example, or in 
combination with one another, could have been used for this purpose. In this case, the 
regression spectrum for TBN was used because that regression spectrum was determined to 
lead to relatively accurate predictions of fluid properties. Other regression spectra may yield 
better results for other data sets. 

[0091] Figure 1 1 shows reduced spectral matrix 1 100. Header row 1 122 lists the points 
in the combined spectrum that have been identified as principal components. Column 1102 
lists the identifiers for the fifteen different lubricants for which combined impedance spectra 
were generated. Note that the fifteen lubricants identified in column 1102 are the same 
fifteen samples identified in column 902 of result matrix 900. Columns 1104, 1106, 1108, 
1 1 10, 1 1 12, 1 1 14, 1 1 16, and 1 1 18 each contain values from the combined plots of the 
resistive and reactive impedance spectra for each of the fifteen lubricants. The identifiers in 
column 1102 are the same as the identifiers in column 902 contained in result matrix 900 
depicted in Figure 9; Figures 9 and 1 1 are based on the same laboratory test of the invention, 
and thus represent a result matrix and reduced spectral matrix respectively based on the same 
lubricating fluid samples. 

[0092] STATISTICAL ANALYSIS OF THE REDUCED SPECTRAL MATRIX 
[0093] As discussed above with respect to block 310 of Figure 3, a number of different 
statistical techniques, some of which use linear relationships and some of which use 
nonlinear relationships, may be used to identify coefficients for use in prediction equations 
(2) above. Linear regression methods with which the present invention may be practiced 
include but are not necessarily limited to Multivariate Least Squares Regression (MLR), also 
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known as Multiple Linear Regression, Inverse Least Squares, or P-Matrix, as well as 
Principal Component Regression (PCR). 

[0094] The present invention takes advantage of the following linear relationship: 


where result matrix R represents a concentration of a property or properties in a fluid, A 
represents a reduced spectral matrix, and coefficient matrix Pisa matrix of loading 
coefficients. The matrix E is an error matrix, also known as the residual matrix, because it 
represents the difference between fluid properties estimated by the expression PA and the 
true value of fluid properties. It will be readily apparent to one skilled in the art that, where 
R and A are known, MLR can be performed to determine P. The matrix P in turn comprises 
coefficients that can be used in prediction equations 124 (1) to obtain predictions of property 
values, such as predicted property values 128, described above with reference to system 100. 
[0101] PCR extends MLR to regress against the scores obtained from PCA as opposed to 
spectral data itself. Recall that the theory of PCA scores matrix S is that S can be used to 
reconstruct the spectral matrix A. Assuming that there is a linear relationship between 
spectral matrix A and concentration C, represented in the present invention by result matrix 
R, it is true that 


It will be readily apparent to one skilled in the art that, where R and S are known, MLR can 
be performed to determine B. The matrix B in turn comprises coefficients that can be used in 
prediction equations (2) to obtain predicted property values 128. 
[0102] Nonlinear regression techniques with which the invention has been practiced 
include pattern recognition analyses, cluster analyses, and neural network analyses. One 
nonlinear regression technique that has been successfully applied in the present invention is 
the Group Method of Data Handling (GMDH) algorithm using the software program 
KnowledgeMiner available on the World Wide Web from Script Software. First, a software 
program such as Mathematica® from Wolfram Research, Inc. of Champaign, Illinois, is used 
to read selected resistive and reactive impedance values for a given set of frequencies from a 
set of files. For example, approximately one frequency at each decade of the frequency 
range could be used. In one instance of practicing the invention impedance values at 10kHz, 
1kHz, 100 Hz, 1Hz, 0.1Hz, 0.01 Hz were selected. 


R = PA + E 


(3) 


R = BS + E 


(4) 
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[0103] The next step is to calculate internal variables that assist in describing a geometric 
shape profile of each impedance spectra. These may consist of the length vector between 
each possible pair of two frequencies that have been collected. The angle, delta, of each 
selected point on the Nyquist representation of the spectra can be included; this angle is 
generally referred to as the dissipation factor, and the tangent of angle delta is calculated by 
dividing resistive impedance by reactive impedance. The tangent of angle delta represents 
the ratio of energy loss versus energy stored at a given frequency. Also, the inverse values of 
the afore-mentioned measured and internal variables are also calculated, allowing for 
variable combinations that include both quotients and products. 
[0104] In order to determine the presence of redundancy in the variables a cluster 
analysis can be preformed. Cluster analysis is a statistical procedure that separates and 
groups a set of data into smaller sets of similar data. If multiple variables are found to be 
statistically similar, one variable from this variable cluster is selected. The resulting reduced 
variable set is then paired with an appropriate physio-chemical parameter to which a 
correlation is to be determined. This data is then modeled using Group Method for Data 
Handling, which is found in commercially available software, i.e., KnowledgeMiner. The 
resulting correlation could be considered a multi-layered neural network composed of 
connection weights that are polynomial, (including linear) functions. 
[0105] Other applications of nonlinear regression techniques to the reduced spectral 
matrix may be apparent to skilled artisans. Further, skilled artisans will understand that the 
aforementioned linear regression techniques will be used to find loading coefficients for 
prediction equations 124 (1) that can in turn be used to obtain predicted property values 128. 
[0106] RUL AND EOL DETERMINATIONS 

[0107] As discussed above with respect to step 312 described in Figure 3, predicted fluid 
properties can be used to estimate either the Remaining Useful Life (RUL) of the fluid, or 
when the fluid will reach its End of Life (EOL). Referring to Figure 1, microcomputer 116 
can be programmed with a variety of logical functions comparing baseline values 122 to the 
results of prediction equations 124 (1), namely property predictions 126, to determine a 
fluid's RUL and/or EOL. 

[0108] As discussed above, each of the baseline values 122 has as attributes the fluid 
property for which the value was recorded, the age of the fluid sample with which the value 
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is associated, and the temperature with which the value is associated. Measurements of RUL 
and EOL may take a variety of forms. For example, RUL may be expressed as a percentage 
of remaining useful life, as a percentage of useful life expended, as the number of miles a 
vehicle can be run before the fluid will reach EOL status, or as the number of hours for 
which an engine can be run before a fluid will reach EOL status. Similarly, EOL may 
represent a value for a fluid property at which the fluid has reached the end of its useful life 
or an age in hours at which the fluid has reached the end of its useful life. 
[0109] In general, microcomputer 1 16 will be programmed with at least one logical 
function that evaluates an expression comparing one or more of the baseline values 122 with 
one or more of the property predictions 126. It should also be understood that the 
determination of RUL and EOL may be based on one or a plurality of fluid properties. 
[0110] In the simplest case, determining RUL or EOL based on one fluid property, 
programming of microcomputer 116 could require iteratively evaluating the expressions (Y > 
vj, Y > v 2 , . . . , Y > v n ) where Y is one of the property predictions 126, each value Vj is one 
of the baseline values 122, and n is the number of fluid samples representing different fluid 
ages for which baseline values 122 have been stored in information library 120. Iterations 
would continue until one of the expressions (Y > vj) evaluated to FALSE, i.e., until the 
program established that prediction of the fluid property value made in real time exceeded a 
threshold value, thus enabling the program to predict the age of the fluid based on the 
threshold value exceeded. 

[0111] In a more general case, programming could iteratively evaluate a compound 
logical expression such as (Yi > vj, \ OR Y 2 > v 2> i OR . . . OR Y k > v kj n ) where Y\ . . . Y k are 
property predictions 126, vi, i . . . Vk, n are baseline values 122, k is the number of fluid 
properties being evaluated, and n is the number of fluid samples representing different fluid 
ages for which baseline values 124 have been stored in information library 120. Thus, if any 
of the fluid properties of interest fail to meet a threshold established by the relevant baseline 
values 122, the expression will evaluate to FALSE, and an RUL determination can be made 
based on the age attribute associated with v\ t - x . . . v k> n . Similarly, microcomputer 1 16 could 
be programmed to make an EOL determination when any given threshold for EOL and/or 
RUL failed to be met. Also, the logical operator OR could be replaced by the logical 
operator AND in the logical expression, requiring that certain thresholds be met for all or at 
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least a subset of property predictions 126 before and EOL and/or RUL determination is 
made. 

[0112] The preceding discussion of the programming of microcomputer 1 16 is meant to 
be illustrative rather than limiting, inasmuch as a skilled artisan would recognize that a 
number of different algorithms could be implemented to make EOL and RUL 
determinations. Any determination of EOL and RUL using property predictions 126 or 
predicted property values 128 would be within the scope of the present invention. 
[0113] Although not shown on Figure 1, microcomputer 1 16 provides output of both 
RUL and EOL determinations. In one embodiment, for example, output of EOL 
determinations is manifested in the activation of a warning light on a vehicle console. In 
another embodiment, an RUL determination is manifested on a digital display embedded in a 
vehicle console. These examples are meant to be illustrative rather than limiting, inasmuch 
as any mode of outputting and displaying EOL and RUL determinations would be consistent 
with the present invention. 
[0114] ALTERNATIVE EMBODIMENTS 

[0115] The preceding disclosures of the invention assume that the spectral matrix is 
created as described in Figure 4. According to the embodiment described in Figure 4, Bode 
plots of resistive and reactive impedance spectra for a plurality of fluid samples are laid 
"head to tail" on a combined plot, and the data from the combined plot is then placed into the 
spectral matrix. Figure 14 describes an alternative embodiment in which data from each of 
the Bode plots of resistive and reactive impedance spectra are processed separately, with 
results then combined for use in predicting fluid properties. 

[0116] Prior to step 1400, steps 400 through 418 will have been performed as described 
above with reference to Figure 4. However, instead of combining Bode plots as described 
with reference to step 420 and then creating a reduced spectral matrix as described with 
reference to step 428, in step 1400 a first spectral matrix is created solely from the resistive 
impedance spectra. Similarly, in step 1402, a second spectral matrix is created solely from 
the reactive impedance spectra. 

[0117] Step 1404 follows the process discussed above with reference to figure 6 to 
perform PCA with respect to the first spectral matrix to create a first reduced spectral matrix. 
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Similarly, step 1406 performs PCA with respect to the second spectral matrix to create a 
second reduced spectral matrix. 

[0118] After step 1406 is complete, the invention can be practiced either by combining 
the first and second reduced spectral matrices and performing a regression on the resultant 
combined reduced spectral matrix, or by performing regressions on each of the first and 
second reduced spectral matrices and combining the results. Accordingly at step 1408 a 
check is done to determine whether regression is to be performed on the first and second 
reduced spectral matrices together or separately. 

[0119] If regression is to be performed on the first and second reduced spectral matrices 
together, they are combined in step 1410. The first and second reduced spectral matrices will 
each have n rows, n being the number of fluid samples being tested, but may have different 
numbers of columns because the number of principal components identified for each of the 
first and second reduced spectral matrices may or may not be the same. The first and second 
reduced spectral matrices are combined simply by placing the data in the two matrices side 
by side, so that the combined reduced spectral matrix has n rows, and has the number of 
columns that is the sum of the number of columns in each of the first and second reduced 
spectral matrices. 

[0120] In step 1412, MLR is performed with respect to the combined reduced spectral 
matrix in the same manner as described above with respect to the reduced spectral matrix. 
Similarly, the scores from the PCAs performed on the first and second spectral matrices 
could have been combined and a PCR or other statistical technique performed on the 
combined scores in the manner described above. 

[0121] If the response to the check done in step 1408 is that regression is to be performed 
on the first and second reduced spectral matrices separately, a first MLR is performed on the 
first spectral matrix in step 1414. In step 1416, a second MLR is performed on the second 
spectral matrix. In step 1418 the results of the first and second MLRs are then used in the 
prediction equations 124 (1) in the manner described above, and the results of the prediction 
equations are then combined to yield final predictions of fluid properties. The results of 
prediction equations might be combined in a number of ways; in one embodiment they are 
averaged. Again, PCR or other statistical techniques also could have been performed in this 
manner. 
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[0122] Embodiments of the invention discussed thus far comprise using at least one Bode 
plot to create a spectral matrix. However, some embodiments of the invention create a 
spectral matrix without using data from a Bode plot, i.e., make a spectral matrix containing 
only Nyquist-derived datum values that are processed using the techniques that have been 
described above with respect to Bode plots. As discussed above with reference to Figure 12, 
there is a considerable amount of data incidental to practicing impedance spectroscopy that is 
not contained in or derived from a Bode plot. For example, as discussed above regarding 
Figure 12, data from Nyquist plots of impedance spectra has also been found to be useful 
when added to a combined Bode plot of resistive and reactive impedance, and then used to 
create a spectral matrix. Some embodiments of the present invention forgo the use of data 
from Bode plots and construct the spectral matrix solely from data such as the data discussed 
with reference to Figure 12. Processing otherwise proceeds as described above with 
reference to Figure 4. 

[0123] The above description is intended to be illustrative and not restrictive. Many 
embodiments and applications other than the examples provided would be apparent to those 
of skill in the art upon reading the above description. The scope of the invention should be 
determined, not with reference to the above description, but should instead be determined 
with reference to the appended claims, along with the full scope of equivalents to which such 
claims are entitled. It is anticipated and intended that future developments will occur in the 
field of determining fluid conditions using impedance spectroscopy and that the disclosed 
systems and methods will be incorporated into such future embodiments. Accordingly, it 
will be understood that the invention is capable of modification and variation and is limited 
only by the following claims. 
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